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ABSTRACT 


A linear combination of modified Bessel functions is defined, 
discussed briefly, and tabulated; namely, 

^(Ax) = I k Ux) K^AXj) - (-l) k+£ I A (A Xl ) ^(Ax) 

This combination was found to recur in the analysis of various heat 
transfer problems and in the analysis of the thermal behavior of liv- 
ing tissue when modeled by cylindrical shells. 
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1 . INTRODUCTION 


Many physics and engineering problems are amenable to solutions 
in terms of Bessel functions. These functions have been extensively 
studied and tabulated. A hitherto undefined linear combination of 
modified Bessel functions was found to recur in the analysis of nu- 
merous engineering problems dealing with diffusion phenomena. Some 
examples in heat transfer are: steady heat flow in thin rods, taper- 
ing fins and thin fins around cylinders , electrical transmission lines 
Laplace transforms applied to flow of heat in cylinders, and the model 
ing of living tissue by the "bio-heat" equation. This combination of 
Bessel functions is defined, discussed, and tabulated in this paper. 



2. ANALYSIS 


Consideration of a steady state energy balance in a heat 
generating material whose heat source is linearly dependent on the 
temperature leads to the following equation: 

V 2 T - A 2 T = 0 (1) 

where T is the temperature of the medium and A 2 is the ratio of the 
strength of the source per unit time, degree and unit volume to the 
conductivity of the medium. In many cases, and for cylindrical 
geometries in particular, Eq. (1) may yield the following equation [1,2]* 

Y* ' + Y' - (A 2 + v 2 /x ) Y = 0 (2) 

Equation (2) is referred to as the modified Bessel equation of order v. 

If V is an integer, the following linear combination of modified 
Bessel functions is found to recur in the solution of Eq. (2): 

^ £ (Ax) e l k (Ax) K^CAXj) - (-l) k+£ I^AXj) yAx) (3a) 

or 

I k (Ax) (-l) k+£ yAx) 

^(Ax) E (3b) 

V Xx i } V Ax i } 

This shorthand definition of the linear combination of modified Bessel 
functions simplifies the analysis of the problems considerably. Further- 
more, it was observed that the function defined in Eqs. (3a) or (3b) 

*Numbers in brackets refer to entries in REFERENCES. 
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exhibits a few interesting characteristics that render the above defi- 
nition even more useful. These characteristics are: 

(a) By definition, the function 'F^(Xx) satisfies the modified 
Bessel equation of order k. 


1 d_ 
x dx 


dx 


X 2 - — W^(Xx) = 0 


(4) 


(b) This function obeys the same differentiation rule as does 
the modified Bessel function of the first kind. 


dx 




kH kil (Xx) + Xx ' 1< (k+l)£ (Ax) 


(5) 


Therefore, it is not a so-called cylinder function. 

(c) For two consecutive indices the following identity is ob- 
tained [3 (p. 375)]. 


T. 


k(k+l) 


(Ax,) = 


Ax. 


( 6 ) 


(d) When the indices are reversed in order, the following ex- 
pression is obtained: 

W = <- i > k+i+1 (7) 

(e) For identical indices the following result is obtained 

Y kk (Xx i )E ° (8) 

The above list of characteristics is not an exhaustive one. 
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When plotted on semi-log coordinates , .the function (A.x) 

appears to behave as two straight lines as shown in Figs. 1 and 2. 
Also, the absolute magnitudes of the slopes of these lines seem to 
be identical. These characteristics suggest the approximation of 
the function by an exponential expression of the kind 


^k(k+l)^ X ) ~ exp U ^ x i ^ + s ^ x 


(9) 


for 


%5 < Ax < X„ 

2 

and 

Ax > X 3 

where X < A x and X > A x . The branches are symmetrical with re- 

2 i o 1 

spect to Ax = Ax t with slopes of S - -0.43 and S - 0.43, respectively. 
For the right-hand branch, i.e.. Ax > X 3 , U is about -0.5AXj , whereas 
for the left-hand branch U is of the order of 0.4Ax 1 . 

Figures 3 and 4 show the function for repeated subscripts 0 and 1. 
A simple computer program was written for the purpose of computing 
the function T^(Ax) = 'F^(X). Results for the first four combinations 
of indices are given in the tables. Values of the arguments, 

X as well as XI, in the tables are for the range 0.01 to 100. 
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TABLES OF THE FUNCTION 

' F kl (X) E V X) k l ( x1) " ( " 1)K+L t l (x1) V x) 
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